Surface rearrangement at complex adsorbate-substrate interfaces. 

E.V. Vakarin and J. P. Badiali 

LECA ENSCP-UPMC, 11 rue P. et M. Curie, 75231 Cedex 05, Pans, France 

On the basis of the information theory approach we propose a novel statistical scheme 
for analyzing the evolution of coupled adsorbate-substrate systems, in which the substrate 
undergoes the adsorbate-induced transformations. A relation between the substrate mor- 
phology and the adsorbate thermodynamic state is established. This allows one to esti- 
mate the surface structure in terms of incomplete experimental information and the one 
concerning the adsorbate thermodynamic response to the structural modifications. 



I. INTRODUCTION 

Adsorbate-induced changes in the substrate 
morphology is a quite common phenomenon in 
nature. This concerns with in-plane and out-of- 
plane surface modifications, such as reconstruction 
or roughening (see [1] for a recent review). 

For simple crystalline surfaces (single crystals) 
there are well-developed Hamiltonian models deal- 
ing theoretically with the adsorbate-induced re- 
structuring effects [2]. Focusing on the roughness, 
the solid-on-solid models [3,4] are appropriate for 
the description of the surface morphology under 
the influence of adsorption. Nevertheless, even 
in this "traditional" domain there are certain un- 
solved problems with the surface structure deter- 
mination [5]. Therefore, one has to work under 
conditions of incomplete information. 

Moreover, many adsorbents do not fall into this 
class because of their intrinsic complexity. The lat- 
ter appears as a consequence of the surface disor- 
der, heterogeneity [6] or polydispersity. Glassy [7], 
granular, amorphous [8], or porous materials, like 
aerogels [9,10], could serve as examples. The mor- 
phology of such materials is conventionally speci- 
fied by a probability distribution (e.g. site ener- 
gies, pore or grain sizes). In many cases the dis- 
tribution itself is unknown. For randomly hetero- 
geneous surfaces the information comes from scat- 
tering experiments or STM images. In the case of 
three-dimensional structures one has to rely upon 
results of indirect measurements (such as nitrogen 
adsorption, mercury porosimetry, thermoporome- 
try, etc). These techniques give a quite limited in- 
formation from which the probability distribution 
should be determined. This problem is beyond the 
scope of the standard statistical mechanics, which 
starts from a microscopic model. 

Actually, the situation is complicated by the fact 
that the surface morphology [11] and/or energet- 
ics [12] change in the course of adsorption. This 
is quite common phenomenon, occurring, for in- 
stance, during the aerogel characterization [9,10] or 
in the course of various insertion processes [13,14]. 
In this respect the problem is essentially as fol- 



lows. One has to describe the surface evolution 
from a poorly characterized initial (e. g. clean 
surface) to a final (adsorbate-induced) state. In 
addition, the adsorbate thermodynamics, coupled 
to the surface morphology, should be calculated. 
The inverse problem, that is, restoring the clean 
surface characteristics from its morphology in the 
covered state, is also important. 

In this paper we make an attempt to deal with 
such problems combining the information theory 
approach [15,16] and the standard statistical ther- 
modynamics. The adsorbate-substrate system can 
be split into subsystems with different levels of de- 
scription. One is a dynamic subsystem (the adsor- 
bate), which evolves according to a Hamiltonian, 
containing unknown parameters (local roughness 
of the surface). The surface is considered as a 
stochastic system. Its morphology (e.g. step con- 
figuration) is specified by a probability distribu- 
tion. This combined scheme looks quite similar 
to the so-called superstatistics [17], that describes 
the statistics in systems with fluctuating intensive 
quantities. The key difference is that in our case 
the probability distribution is unknown. It has to 
be inferred from a limited number of constraints. 
Recently such a maximum-entropy procedure has 
been applied to a description of velocity fluctua- 
tions in turbulent fluids [18] and the strain fluctu- 
ations in heterogeneous intercalation systems [19]. 
In our case the adsorbate can influence the sur- 
face roughness and their coupling determines both 
the thermodynamic behavior and the shape of the 
probability distribution. It is essential that one 
necessarily deals with two entropic impacts [20]: 
the thermodynamic entropy (due to the adsorbate) 
and the information entropy (due to the proba- 
bility distribution). Constructing from these two 
terms a suitable entropy measure, we investigate 
the evolution of the surface morphology and the 
adsorbate thermodynamics. The main purpose of 
this paper is to find a relation between the sur- 
face structure and the adsorbate thermodynamic 
response to these structural modifications. This 
would allow one to gain an information on the sur- 
face following the adsorbate behavior. 
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II. CLEAN SURFACE 

In this section we demonstrate how the maxi- 
mum entropy formalism can be applied to inferring 
a distribution based on a limited information. For 
concreteness we consider a solid-like semi-infinite 
substrate (see Figure 1). The surface (substrate) 
morphology is represented by a collection of steps 
{hi} growing on discrete sites i = l..N s , with 
— oo < hi < oo, describing a deviation from the 
flat geometry. A sketch of the surface profile is 
given in Figure 1. 

It is assumed that, because of its complex na- 
ture, the surface cannot be described determinis- 
tically, i.e. in terms of a Hamiltonian governing 
the substrate evolution. Thus we have to accept 
a probabilistic description, focusing on a distribu- 
tion P({hi}) of the step configurations. However, 
the distribution is unknown. The only available in- 
formation concerns with a set of experimental tests 
determining [21], for instance, the average rough- 
ness 



H = (dhi)P({hi})^2h^ 



(1) 



For simplicity the steps are taken as spa- 
tially non-correlated P({hi}) — YliP(hi). Then 
^2ihf = N s h 2 , where N s is the number of steps. 
Therefore, P(h k ) is a probability to find a step of 
magnitude ft at a randomly chosen surface site k. 

Having such a restricted information on the sub- 
strate, one can infer its statistics within the in- 
formation theory approach [15,16]. The inference 
procedure is based on a given entropy measure Si, 
which is taken here in the Shannon form 



Si = -J2 dh k P(h k )\nP(h k ) 



(2) 



This entropy is a measure of our uncertainty on 
the system, in the sense that Si increases with 
broadening P(h k ). The least biased normalized 
distribution, compatible with the constraint (1) is 
obtained by extremizing the Lagrangian 



A = Sj-i/^ Jdh k P(h k )-l 

k L 

]T J dh k P{h k )h 2 k - H Cj 



(3) 



A 



This leads to a gaussian distribution 
Po(hi 



The Lagrange multiplier A can be determined 
from the constraint (1) Ao = 1/(2Hq). We do not 
discuss here how precise is this description. Hav- 
ing more information, e.g. (h 4 ), we can follow the 
same scheme in order to obtain a refined distribu- 
tion. 



III. SURFACE-ADSORBATE COUPLING 

Being exposed to a fluid, the surface can adsorb 
the fluid particles on the sites (see Fig. 1). Having 
information [21] on the roughness H in the pres- 
ence of adsorption, one can infer the step distribu- 
tion maximizing Si. Obviously, the distribution is 
also gaussian 



P(hi) 



J dhie 



-\h 2 



(5) 



/ dhie- Xoh i 



(4) 



but with a different width A = 1/(2H). There- 
fore, if H is close to H n , then one would conclude 
that P(hi) = Po{hi) and the adsorbate-induced ef- 
fects are negligible. This can be the case if the 
surface is quenched. As is discussed above, the 
surface may change its morphology even if the 
average roughness remains unchanged. In other 
words, we do not have new a priori information if 
H = Hq. Then, starting with the same entropy 
measure Si , we arrive at a trivial (but misleading) 
result: P(hi) = Po(hi). This defect can be avoided 
in (at least) two ways. One might need to compare 
higher moments (for instance, the average of hj). 
This would require additional tests [21]. Alterna- 
tively, we may refine the entropy measure, taking 
into account the adsorbate-induced effects. This 
requires to make a link between the surface struc- 
ture and the adsorbate properties. 

The adsorbate is considered as a dynamic sub- 
system, whose coupling to the surface is given by 
a Hamiltonian 7i({ti}\{h k }). We assume that the 
surface evolves in time much slower than the adsor- 
bate, that reaches the equilibrium while the surface 
is practically unchanged. Tracing over the set of 
occupation numbers {U} at a given step configu- 
ration we can calculate the partition function 

Z(T\{h k })=Y,e- mAU]l{hk}) 
{U} 

where (3 = is the inverse temperature, and 

r is a thermodynamic variable (e. g., chemical po- 
tential or coverage, see below) . Then the adsorbate 
thermodynamics, for instance, the conditional en- 
tropy ST(r\{h k }), can also be calculated. 

In what follows this quantity will play an essen- 
tial role. Therefore, it is important to clarify its 
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meaning. As a function of r the entropy is the 
conventional thermodynamic state function. It is 
a measure of the adsorbate disorder for a given 
surface structure {h k }- On the other hand, if r 
is fixed (e. g. fixed adsorbate pressure), then St 
as a function of {h k } should be considered as an 
effective potential, selecting a preferable step con- 
figuration the one which maximizes St- 

Since we are dealing with adsorbate-induced ef- 
fects, then ST(r\{hk}) should be relevant to the 
surface structure evolution. Therefore we refine 
the entropy measure 



Sj + kJ (dh k )P({h k })S T (r\{h k }) 



(6) 



which takes into account the fact that some part 
of our uncertainty on the surface state is due to 
the adsorbate thermodynamic state 5r(r|{/ifc}). 
In what follows we assume that the steps do not 
correlate through the adsorbate Sr^K^fc}) = 
X^fc ST(i~\hk). Imlicitly this means a restriction on 
the Hamiltonian 7i (see below). 

The parameter k is a measure of the adsorbate 
influence on the step distribution. In other words, 
k is a degree of the surface quenching, such that 
S — ► Sj as k — » and we return to the clean 
surface problem. Determination of k is a separate 
problem that might require additional tests. For 
instance, comparing the surface step creation en- 
ergy e and the adsorbate binding energy excess 7, 
wc may estimate n = j/(e + 7). However, in this 
study k is taken as a parameter. 

Therefore, we consider the Lagrangian 



A' = A + K 



dh k P(h k )S T (T\h k ) (7) 



that must be extremized. This procedure leads to 

-X K h^+KS T (T\hi) 



P K (hi) = 



J dhie-^ h i+ KS ^ T ^ 



(8) 



The multiplier A K should be determined from the 
constraint on the roughness H. The term A K /i| 
partially keeps traces of the surface preparation (a 
given roughness H). The second term KST(r\hi) 
corresponds to the adsorbate-induced effects. If St 
is an analytic function, then we may expand [22] 

S T {r\h l )=Y,A n {r)h- 

n 

Therefore, even without resorting to any specific 
model, it is clear that P K (hi) is potentially more in- 
formative than P{hi) as the former contains higher 
order terms (in hi). As is noted above, in or- 
der to obtain a similar result staring from Si, one 



needs additional information. Therefore, refining 
the entropy measure, we restore (at least partially) 
this information. In particular, we can calculate a 
global change of the surface structure in compar- 
ison to the clean surface case. This can be done 
using the symmetric Kullback-Leibler information 
measure 



K = Y,jdh l [P K {hi) - P {hi)]\n 



P K {hi 



Po(hi) 



(9) 



which tells us how much these two distributions 
are different. K is found to be 

K = -(\ K - \ )(H - H ) (10) 
+K^2[(Sr(T\hi)) K - {S T (r\hi))o] 

i 

where the angular brackets (...) and (...)« de- 
note the averages taken with the distributions 
Po(hi) and P K (hi), respectively. Therefore, even 
if H = Hq, the global variation of the surface 
structure, induced by adsorption is proportional to 
the change of the adsorbate thermodynamic state. 
The state (...)o does not necessarily correspond to 
the clean surface, it could be any suitable reference 
state. In the same way one can estimate the infor- 
mation gain due to the refinement of the distribu- 
tion (replacement P(hi) by P K (hi)). This result 
is quite general, it does not depend on the micro- 
scopic details underlying ST{r\hi). Nevertheless, 
in order to study the distribution itself we need to 
consider a microscopic model. 



IV. ILLUSTRATION 

In order to demonstrate what kind of local 
changes are induced by adsorption, as an illustra- 
tion we choose one of the simplest adsorption mod- 
els - the non-interacting lattice gas with a height- 
dependent binding energy [23,24]) 

H({t i }\{h k }) = -J2(» + e(h i ))t i (11) 



We are aware that the model itself is quite " rough" 
in application to real adsorbates. It neglects the 
lateral interactions and, therefore, is inadequate 
for the description of the adsorbate critical be- 
havior (such as liquid-gas or order-disorder transi- 
tions). Nevertheless above the critical temperature 
this model is qualitatively correct. In addition it 
is exactly solvable such that one can be sure that 
our implementations arc not artifacts of approxi- 
mations (e.g. mean field). Moreover, we do not 
take into account that the roughness increases the 
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surface area, allowing for more than a monolayer 
coverage [24]. Thus t, = or 1 independently of 
hi. 

The entropy of this simplified model is given by 
(the role of r is played by the chemical potential 

S T {ii\hi) = -6{hi) H6(hi)] - [1 - 0(hi)] ln[l - 0(hi)] 

(12) 

Here the coverage is Langmuirian 

Le 0e( hi ) 



e(hi) = 



1 + LeP< h >"> 



(13) 



The adsorbate activity L is related to the chemi- 
cal potential [i via L = exp(/3/i), and e(hi) is the 
height-dependent adsorbate binding energy. 



e(h) =e + 7/i? 



(14) 



A step can be viewed as a loosely coordinated site 
with the binding energy larger [23,24] than that 
for the flat surface. Here 7 is the binding energy 
excess due to the local deviation from the flat ge- 
ometry. The binding to the flat surface eo is an 
irrelevant constant that can be absorbed into the 
definition of L. 

In Figure 2 Sri^hi) is analyzed as a function 
of hi and dimensionless activity m = + eo). 
The entropy (as a function of m) behaves in the 
usual way, exhibiting a single maximum at some 
m* which depends on the magnitude of hi and 
7. At relatively high adsorbate activities (m > 0) 
Sriplhi) is also single-peaked function of hi, 
with the maximum located at hi = 0. This sug- 
gests a tendency towards random deviations from 
hi = 0. This maximum transforms into a minimum 
when the activity changes its sign (m < 0). Then 
the entropy develops two symmetric peaks at some 
\h*\ that depends on the magnitude of m and 7. 
This means that the position hi = becomes un- 
favorable and the surface would tend to exhibit a 
random zig-zagging with the steps distributed near 
±h*. Therefore, considering Sri^hi) as an effec- 
tive potential for the step variables hi , one expects 
qualitative changes in the surface morphology with 
changing adsorbate activity m. 

In order to analyze (qualitatively) what this im- 
plies for the distribution P K (hi), we take A K as 
a parameter and set the normalization constant 
(the denominator in (8)) equal 1. Then we can- 
not calculate averages, but these simplifications do 
not change the shape of the distribution. Such 
a reduced P K (h) is plotted schematically in Fig- 
ure 3. It is seen that with increasing k the dis- 
tribution changes from a gaussian (k = 0) to a 



tri-modal. This results from an interplay between 
the adsorbate-induced effects and the tendency to 
keep a given roughness. The curves coincide in the 
\&rge-h region because the distribution is not nor- 
malized. For quantitative purposes this defect can 
be removed by numerical integration in eq. (8). It 
is clear (see Fig. 1) that the adsorbed species oc- 
cupy the space of the corrugated surface layer [23] . 
Therefore the adsorbate is confined to a non-planar 
slit, whose average thickness T is determined by 
the surface roughness H (see Figure 1). Moreover, 
if hi gives a length scale, then the surface roughness 
H (the average of hj) is at the same time the 
surface area excess in comparison to a flat surface. 
Thus P K (hi) can be considered as a surface area 
distribution. Within this analogy we make a qual- 
itative comparison to Monte Carlo results [11] for 
the surface area fraction in pores with rough walls 
(see the insets in Figure. 2). The right inset ex- 
hibits a peaked (nearly gaussian) distribution that 
corresponds to an empty pore. The left inset de- 
picts the adsorbate-induced surface structure. It 
is seen that our approach describes quite well this 
physics (recall that n is a measure of adsorbate- 
induced effects). We did not try to fit, because 
the models are not strictly identical. It is remark- 
able that our scheme is capable of inferring the 
surface structure even without entering the micro- 
scopic details concerning the surface. Moreover the 
simulation results [11] are not much sensitive to 
the type of adsorbate- wall interaction (hydrophilic 
or hydrophobic, see Fig. 5 in [11]). This suggests 
that the adsorbate-induced surface rearrangement 
is mainly an entropic effect. This is coherent with 
our implications. 



V. CONCLUSION 



The evolution of coupled adsorbate-substrate 
system is investigated in terms of a new approach 
combining the information theory and the stan- 
dard statistical mechanics. The approach is suit- 
able for complex systems, where, for a number of 
reasons, not all the microscopic details (e. g. a 
Hamiltonian) are available. 

The approach establishes a clear correlation be- 
tween the substrate morphology (e. g. roughness) 
and the adsorbate thermodynamic state. This 
implies a statistical scheme for surface character- 
ization in the course of the adsorbate-substrate 
coupling, and allows one to estimate the surface 
morphology based on the information concerning 
the adsorbate thermodynamic response. One can 
also predict the clean surface properties based on 
its morphology under adsorption and vice versa. 
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Thus, our scheme is, in some sense, complimen- 
tary to the recent (mainly simulational) surface 
characterization techniques [6,11,12]. It can also 
be applied to analyzing the experimental data on 
porous adsorbents [9,10]. 

We would like to underline that the simple model 
considered here serves just as an illustration. To 
be more realistic in further applications we have to 
consider the adsorbate-adsorbate interactions and 
a correlation between the steps. Also, the exper- 
imentally available information (H or H) is sup- 
posed to be quite crude. In many cases of real sur- 
faces a more detailed statistics can be extracted 
(for instance, from STM images), for example, - 
pair correlation Hij — ((hi — hj) 2 ). Nevertheless, 
for three-dimensional matrices (like aerogels) in 
the majority of cases only indirect measurements 
are available. All these features can be easily in- 
corporated into our approach. 

Apart from the practical applications, our ap- 
proach seems promising in analyzing a class of 
coupled dynamic-stochastic systems, where the 
stochastic subsystem is not perfectly quenched. 
Particularly interesting is to study the adsorbate 
critical behavior under conditions when a porous 
matrix responds to the adsorbate thermodynamics 
(e.g. helium absorption in aerogels). 
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FIG. 1. A sketch of the adsorbate-substrate interface. The black dots mark the adsorption sites. The surface 
steps hi are defined for each site. The solid line is a reference from which all hi are counted. 



FIG. 2. The conditional entropy Srinlhi) as a function of hi and dimensionless activity m = (3([a + eo), /?7 = 1 



FIG. 3. A non-gaussian step probability distribution, resulting from eq. (8) (A K = 0.1, L — 0.1, f3f = 2). The 
insets demonstrate the Monte Carlo results [11] for the surface area fraction in pores with rough walls. 
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